Mon. Not. R. Astron. Soc. 000, [|-?? (2002) Printed 1 February 2008 (MN 1*1^ style file vl.4) 



Type la supernovae and the formation history of 
early-type galaxies 

Ignacio Ferreras and Joseph Silk* 

Physics Dept. Denys Wilkinson Building, 1 Keble Road, Oxford 0X1 3RH, United Kingdom 



Submitted March 18, 2002; Revised June 11, 2002; Accepted July 17, 2002 



(N 
O 
O 



0^ 
(N 

> 
(N 
cn 

o 

(N 
O 

o 



ABSTRACT 

Using the standard prescription for the rates of supernovae type II and type la, we 
compare the predictions of a simple model of star formation in galaxies with the 
observed radial gradients of abundance ratios in a sample of early-type galaxies to infer 
the relative contribution of each type of supernova. The data suggests a correlation 
between the fractional contribution of Type la to the chemical enrichment of the stellar 
populations (1 — £,) and central velocity dispersion of order 1 — ^ ~ —0.16 log ctq + 0.40, 
so that the type la contribution in stars ranges from a negligible amount in massive 
(ao ~ 300 km s~-^) galaxies up to 10% in low-mass 100 km s~^) elliptical galaxies. 
Our model is parametrized by a star formation timescale (tsp) which controls the 
duration of the starburst. A correlation with galaxy radius as a power law (igp (x r^) 
translates into a radial gradient of the abundance ratios. The data implies a wide range 
of formation scenarios for a simple model that fixes the luminosity profile, ranging from 
inside-out (/? = 2), to outside-in formation (/3 = —1), as is consistent with numerical 
simulations of elliptical galaxy formation. An alternative scenario that links tsF to the 
dynamical timescale favours inside-out formation over a smaller range 0.4 < (3 < 0.6. 
In both cases, massive galaxies are predicted to have undergone a more extended period 
of star formation in the outer regions with respect to their low-mass counterparts. 
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1 INTRODUCTION 

Type la supernovae (SNIa) describe stellar explosions whose 
spectra show lines of elements or intermediate mass such as 
silicon, and of the iron group, but no hydrogen lines. They 
currently occupy a very special place in cosmology since they 
can be used as standard candles. Even though their lumi- 
nosities span an order of magnitude, empirical correlations 
between their absolute magnitude and the shape of their 
light curves can be used in order to determine distances 
(Phillips 1993) in an analogous way to classical Cepheids. 
Furthermore, the absolute luminosity of SNIa is orders of 
magnitude brighter than variable stars, enabling us to ob- 
serve them at high redshift. This is a tool currently being 
used to constrain the cosmological parameters (e.g. Perlmut- 
ter et al. 1999; Riess et al. 1998). However, SNIa also play 
a very important role in the process of galactic chemical en- 
richment since they contribute a large amount of iron to the 
interstellar medium (Thielemann, Nomoto & Yokoi 1986). 
In fact Type la supernovae might even be the main produc- 
ers of iron in the Universe (Ishimaru & Arimoto 1997; how- 
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ever see Gibson, Loewenstein & Mushotzky 1997). The ob- 
served properties of SNIa's suggest a binary model in which 
at least one of the stars is a white dwarf that reaches the 
Chandrasekhar limit (~ I.4A/0) by accretion or by merging 
with another white dwarf. The timescale for the explosion 
is thereby limited by the lifetimes of stars which end up as 
white dwarfs, i.e. with masses M ^ 8Mq. This implies that 
Type la's can be observed long after star formation sub- 
sided. Indeed, all supernovae observed in early-type galaxies 
— which feature no ongoing star formation — are Type la, 
whereas late-type galaxies display a mixture of Type la and 
Type II supernovae (Cappellaro et al. 1997) 

On the other hand, Type II supernovae (SNII) show 
hydrogen lines in their spectra and arise from the core col- 
lapse of a single, massive (M ^ 8Mq) star at the end of 
its lifetime, which occurs between 1 — 50 Myrs after the 
core hydrogen burning phase started. Hence, the timescales 
of either type of supernovae are remarkably different. Fur- 
thermore, the yields of chemical elements are also in sharp 
contrast, since SNIa produce much more iron than SNIIs, so 
that stars born during the first phases of star formation — 
when the contribution from SNIa to the interstellar medium 
was negligible — display an enhancement of a elements over 
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iron, with respect to the younger generations of stars, such 
as the Sun, which are born in an environment polluted by 
both types of supernovae. 

Observations of radial gradients in the colours (Franx, 
lUingworth & Heckman 1989; Peletier et al. 1990; J0rgensen, 
Franx & Kjaegaard 1995) and line indices (Gonzalez 1993; 
Davies, Sadler & Peletier 1993; Peletier et al. 1999) in el- 
liptical galaxies show gradients, being mostly redder and 
more metal rich at their centres, although some early-type 
galaxies display blue cores (Menanteau, Abraham & Ellis 
2001). Broadband photometry has been a technique repeat- 
edly used to infer the star formation history of galaxies 
notwithstanding the age-metallicity degeneracy (Worthey 
1994), which prevents us from getting a well-defined pic- 
ture of galaxy formation using colours alone. A combined 
analysis of line indices seem to break the degeneracy since 
their age and metallicity dependence can be rather differ- 
ent (Kuntschner 2000; Trager et al. 2000a, 2000b). Abun- 
dance ratios represent an alternative observable since the 
timescales for SNIa and SNII are remarkably different. Gi- 
ant early-type galaxies feature an overabundance of a ele- 
ments over iron (Peletier 1989; Worthey, Faber & Gonzalez 
1992; Trager et al. 2000a; Kuntschner 2000), indicative of a 
short duration of the star formation stage, so that mostly 
SNII contribute to the metallicity of the stellar component. 

This has been used in several models of galaxy for- 
mation in order to constrain the star formation history of 
galaxies. Matteucci (1994) analysed the observed [Mg/Fe] 
ratios in ellipticals to constrain the star formation history 
to timescales shorter than 0.1 Gyr. Furthermore, the trend 
of [Mg/Fe] with galaxy mass was used to imply either an 
increasing star formation efficiency with galaxy mass, or a 
top-heavy initial mass function, so that more massive stars 
were formed in ellipticals compared to a more quiescent en- 
vironment such as in a disk galaxy. A detailed analysis of 
abundance ratios in metal-poor stars is a valuable tool for 
determining the star formation history in the solar neigh- 
bourhood (Gratton et al. 2000), and this can be extended 
to stellar populations in globular clusters. In a recent anal- 
ysis of abundance ratios in a sample of six red giant stars 
in the to Cen cluster, Pancino et al. (2002) found evidence 
for the contribution of SNIa to the composition of younger, 
more metal rich red giants. 

In this paper we explore the radial gradients observed 
in abundance ratios such as [Mg/Fe] in elliptical galaxies 
as a function of the star formation timescale. The relative 
contribution from SNIa and SNII can be compared with ob- 
servations to infer the formation process of the stellar com- 
ponents, enabling us to connect the star formation history 
and the dynamical history of early-type galaxies. In §2 and 
§3 we describe the model used to predict the rates of either 
type of supernovae and their contribution to chemical en- 
richment. §4 describes our model predictions and compares 
them to observed data. Finally, in §5 we discuss the implica- 
tions of the comparison between our simple model and the 
observed data. 



2 RATES OF TYPE lA AND TYPE II 
SUPERNOVAE 

The progenitors of type la supernovae are still a matter of 
debate. The most common candidates are either the dou- 
ble degenerate model (Iben & Tutukov 1984) in which two 
C-O white dwarfs merge, reaching the Chandrasekhar mass 
and exploding by C deffagration, and the single degenerate 
model (Whelan & Iben 1973) comprising a close binary sys- 
tem made of a nondegenerate star and a C-O white dwarf 
which accretes material from the other component, reaching 
the Chandrasekhar mass and triggering C deflagration. The 
negative result in the search for close binary systems made 
up of two massive white dwarfs in a sample of 54 white 
dwarfs (Bragaglia et al. 1990) hints at the single degenerate 
model as the most plausible one, which is the one we will 
assume henceforth. The ESO SNIa progenitor survey, which 
will use UVES/VLT (Koester et al. 2001) on a sample of 
1,500 white dwarfs will clarify this point. 

In order to estimate the rates of type la supernovae, we 
follow the prescription of Greggio & Renzini (1983), recently 
reviewed by Matteucci & Recchi (2001). The rate can be 
written as a convolution of the initial mass function (IMF) 
over the mass range that can generate a type la supernova. 
The lower mass limit is Msm = 3Mq in order to generate 
a binary with a white dwarf which will reach the Chan- 
drasekhar limit. The upper mass limit is Mbm = 16Mq so 
that neither binary undergoes core collapse. The rate is thus; 

Rla{t) ^ A / dm(l>{m) I f{^')lp{t - Tm2), (1) 

where rrit = max(3MQ, m(t)) and m(t) is the mass of a star 
whose lifetime is t, whereas Tm^ is the lifetime of the non- 
degenerate companion, with mass m2. f(p) is the fraction 
of binaries with a mass fraction ^ = M2/MB, where M2 
and Mb are the masses of the secondary star and the bi- 
nary system, respectively. The range of integration goes from 
Urn = max(m2/m, 1 — 8M0/m), to ^max = 0.5. The analysis 
of Tutukov & Yungelson (1980) on a sample of about 1000 
spectroscopic binary stars suggests that mass ratios close to 
^ — 1/2 are preferred, so that the normalized distribution 
function of binaries can be written: 

/(m) = 2^+^(1 + 7)m^ (2) 

as suggested by Greggio & Renzini (1983), and we adopt the 
value of 7 = 2. The normalization constant A is constrained 
by the ratio between type la and type II supernovae in our 
Galaxy that best fits the observed solar abundances. We use 
the result of Nomoto, Iwamoto & Kishimoto (1997), namely 
Ria/Rii = 0.12 to find A = 0.05. 

The rate of type II supernovae is also given as an inte- 
gral of the star formation convolved with the IMF, in this 
case over the mass range expected for the precursors of sin- 
gle component supernovae, i.e. 8 < M/Mq < 50. Here we 
adopt 50Mq as the upper mass cutoff. We use throughout 
the paper either a Salpeter (1955) or a Scalo (1986) IMF. 

/•50Mq 

Rii{t)= dm<t){m)ij{t-r,^), (3) 

JsMq 

where we have neglected the contribution from binaries to 
the IMF over this mass range. 

Figure 1 shows the rates of either type of supernovae for 
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TABLE 1: Supernovae yields (in Mq) 




Type 


Mg 


a 


Fe 


[a/Fe] 


SNIa 


0.016 


0.438 


0.742 


-1.33 


SNII(Salpeter) 


0.122 


3.177 


0.094 


+0.43 


SNII(Scalo) 


0.102 


2.663 


0.105 


+0.30 



t/Gyr 



Figure 1. Relative contribution (top) and rates (bottom) of 
type la and type II supernovae for a generic star formation his- 
tory described by a gaussian function peaked at 3 Gyr with a star 
formation timescale of tgpO.l (left) and 1 Gyr, respectively. 



a generic star formation history given by a gaussian profile 
with a width of tsF =0.1 Gyr {left panels) or 1 Gyr, respec- 
tively. The shape of the curve plotting the rate of type lis 
is nearly identical to the star formation rate given the very 
short lifetimes of massive stars that undergo core collapse. 
Hence, the rate of this type of supernovae is negligible a few 
Myr after the star formation has ceased. On the other hand, 
the longer timescales of the nondegenerate companions in 
type las imply that these supernovae start contributing sig- 
nificantly to the chemical enrichment of the ISM a few hun- 
dred Myr after the main burst of star formation. This effect 
is more extreme in sharply peaked bursts of star formation 
although its contribution to chemical enrichment is only im- 
portant for the interstellar medium and not so much for the 
metallicities of the stellar populations unless later bursts of 
star formation allow the ejected metals from SNIas to be 
locked into subsequent generations of stars. 



3 CHEMICAL ENRICHMENT 

The supernova rates allow us to infer the amount of metals 
ejected to the interstellar medium. We target magnesium 
and iron as the main elements which should be observed to 
discern the contribution between SNIa and SNII. Magne- 
sium is one of the a elements, which are those that can be 
obtained from the reaction of C and O with helium nuclei. 
These elements are: Ne, Mg, Si, S. The solar abundance of 
these elements is slightly enhanced with respect to the mean 
(given by some average metallicity Z). Another family of el- 
ements — the so-called iron-peak elements, which include Fe 



and Co — are less abundant with respect to Z. Notice that 
the amount of iron-peak elements contributes only a small 
fraction to the total metallicity (8% at solar abundance), 
so that changing the abundance of Fe by a large amount 
does not affect significantly the average metallicity. Hence, 
it would be more correct to think of a depression in Fe-peak 
elements in solar abundances, rather than an enhancement 
of Q elements at the centres of elliptical galaxies (Trager et 
al. 2000b). 

The yields from SNIa are taken from the single de- 
generate model W7 of Thielemann et al. (1986), which as- 
sumes a progenitor made up of a IMq, C+O white dwarf 
with equal parts of ^^C and ^^O, with an accretion rate 
of 4 X 10"** Mq yr"\ The yields from SNII are taken from 
Thielemann, Nomoto & Hashimoto (1996) for a range of 
progenitor masses. The stellar initial mass function (IMF) 
is used in order to take a weighted-average of the yields with 
respect to mass. The yields used in this paper are shown in 
Table 1. Notice that a very significant amount of iron is 
produced in SNIa with respect to the a elements. The mass 
contained in 0-, Ne- and C-burning shells is too small com- 
pared with the mass in the Si-burning zone. This implies 
SNIa ejecta are dominated by the products of complete and 
incomplete Si-burning, i.e. a higher iron yield compared to 
the ejecta from core-collapse supernovae. However, this is 
still debatable because the yields may depend rather sensi- 
tively on the initial metallicity. For instance, assuming an 
average metallicity of 0.5^© over the star formation history 
of the galaxy will lead to a total reduction of 25% in the 
amount of Fe ejected from SNIa (Thielemann et al. 1986). 

In this paper we do not consider the contribution of low 
and intermediate- mass stars (M < 8Mq) to chemical enrich- 
ment. We focus on the contribution from either type of su- 
pernova to the star formation history in early-type galaxies 
through the analysis of Mg and Fe. These elements are not 
produced in stars which do not undergo core-collapse, ex- 
cept for very small traces of Mg in intermediate mass stars 
(Marigo, Bressan & Chiosi 1996), which can be neglected 
for our analysis. However, one could argue that the metal- 
licity dependence of the yields from both supernova types 
will introduce a dependence on the contribution to chemical 
enrichment from lower mass stars. We do not consider this 
dependence in our models but we emphasize this point as a 
possible factor to improve upon as more developed theoret- 
ical models of supernova explosions become available. 

Another important factor that needs further work is the 
"starting" abundance ratios one would expect from Popu- 
lation HI stars, the very first stars formed. They are as- 
sumed to have a very different IMF, so that their masses 
are much higher than subsequent generations of stars. The 
study of the yields from this population is still under way, 
but some preliminary modelling of high-mass helium cores 
(64 — 133Mq) suggests the yields present solar abundances 
for nuclei with even nuclear charge (Si, S, Ar) and signif- 
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icant defRciencies in odd-charge nuclei (Na, Al, P) (Heger 
& Woosley 2002). If the pre-enrichment from Population III 
stars is important, one would have to consider three different 
contributors to the abundance of Mg and Fe. 

Assuming only two contributors, namely SNIa and 
SNII, we can find the correlation between [Mg/Fe] and the 
relative contribution of each type of supernova, quantified 
by the parameter 



Rii 



Rii + Ria 



(4) 



where Ria and Rn are the rates of SNIa and SNII, respec- 
tively. Figure 2 shows [Mg/Fe] as a function of (1 — for two 
different IMFs: Salpeter (<;/>(m) oc m"^'^'^; Salpeter 1955) 
and Scalo (0(m) oc m~^'^^; Scalo 1986). The data points 
are from the detailed analysis of early-type galaxies in the 
Fornax cluster by Kuntschner (2000), and from the sample 
of Trager et al. (2000a), evaluated at two radial positions, 
r = -Re/8 and Re/2. The translation from line indices to 
metal abundances is done by matching the data against a 
grid of simple stellar populations over a range of ages and 
metallicities. The data points are plotted as a function of 
central velocity dispersion (ctq, top axis), which has been 
chosen so that the model prediction — plotted against (1— C) 
— gives a good fit to the data. This implies a correlation of 
order: 



-0.16 log (JO +0.43. 



(5) 



The vertical line at 1 — ^ = corresponds to the predic- 
tion for a model in which no SNIa contribute to the enrich- 
ment of the stellar populations. Hence, the contribution of 
SNIa range from a negligible amount for massive galaxies 
(cro ~ 300 km s"^) up to 10% for low- mass elliptical galax- 
ies ((To ~ 100 km s~^). The timescales associated to the 
triggering of SNIa imply star formation occured in a very 
short burst in massive ellipticals, whereas low-mass galaxies 
have a slightly more extended period of star formation. This 
can be accomodated in the framework of hierarchical clus- 
tering if we assume that the latest merging stages of massive 
ellipticals do not give rise to star formation but, rather, in- 
volve the merging of the stellar and gas component. Solar 
abundances are obtained for a contribution of SNIa to the 
total rate (1 — ^) of order 10%, which should not come as a 
surprise since this is the constraint used in order to compute 
the proportionaly constant A described in §2. 



4 ABUNDANCE GRADIENTS 

Observations of colours (Franx et al 1989; Peletier et al. 
1990; J0rgensen et al. 1995) and line indices (Gonzalez 1993; 
Davies et al. 1993; Peletier et al. 1999) in elliptical galaxies 
show radial gradients, being redder and richer in metals at 
their centres. The observed spectral indices display a metal- 
licity gradient of A[Fe/H]/A log = -0.2 ± 0.1 (Davies et 
al. 1993), consistent with the observed colour gradients. A 
metallicity gradient can be directly linked to the standard 
picture of galaxy formation by the dissipative collapse of a 
gas cloud. Carlberg (1984) analyzed such models with an 
N-body code, and predicted abundance gradients around 
—0.5 in massive galaxies, flattening towards zero slope in 
lower mass galaxies. Recent merging events in massive early- 



log c7j,(km s"') 
2,60 2,40 3,20 3,00 1,80 1,60 




Figure 2. Contribution from SNIa to chemical enrichment. (1 — 5) 
denotes the fractional contribution to the total supernova rate 
from SNIa, so that 1 — ^ = represents a pure contribution from 
corc-coUapse supernovac. The data points are from Trager et al. 
(2000a; TOO), evaluated at Re/8 {GUcd S(3uares) and Re/2 (hollow 
circles), and from Kuntschner (2000; KOO). The observed data are 
plotted against central velocity dispersion (top axis). The model 
prediction is compatible with the data for a linear correlation 
between central velocity dispersion and SNIa contribution such 
as: (1 — 5) = —0.16 log (TO -I- 0.43. The steeper slope of the Scalo 
IMF at the bright end implies a lower contribution from high 
mass stars (A/ ^ 20Mq), which results in a lower enhancement. 



type galaxies will dilute any radial signature from an ear- 
lier monolithic collapse process, thereby flattening the slope 
(White 1980; Mihos & Hernquist 1994). Hence, abundance 
gradients allow us to quantify the degree to which hierarchi- 
cal clustering has contributed to the formation of galaxies. 

We assume a simple model with spherical symmetry and 
a radially-dependent star formation efficiency. The standard 
model of star formation assumes a power-law dependence 
between the star formation rate (t/j) and the gas density 
(pg; Schmidt 1963), namely: 



lll{t) = CeftPg{t) (X exp{-t/tsF), 



(6) 



where the constant of proportionality (Ces) is the star for- 
mation efficiency, and the exponential function is the solu- 
tion for a linear Schmidt law (n = 1), in which the star for- 
mation timescale tsF can be written as a function of the ef- 
ficiency timescale {res oc 1/Ccs) and the gas infall timescale 
(r/) as: 



1 

tsF 



— + — 

TcB Tf 



(7) 



Star formation is known to take place in clouds of molecular 
hydrogen. Cloud-cloud collisions contribute to the collapse, 
cooling and subsequent formation of stars. In this simple 
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Figure 3. [Mg/Fe] profile as a function of the star formation 
timescale (tsp) wliicli parametrizes the generic star formation 
history used in this paper. Several curves are given for different 
ages of the galaxy: 0.5, 1, 5 and 11 Gyr. The numbers in paren- 
thesis give the slope of the profile (A [Mg/Fe] /A log tgp) for an 
age of 11 Gyr. 



scenario, we would expect the star formation efficiency to 
be proportional to the frequency of cloud-cloud collisions, 
which should scale with the number density of clouds. Hence, 
the gas density profile traces the star formation efficiency, 
so that at the centre the star formation timescale should be 
shorter than in the outer regions. In this scenario tsF is a 
monotonically increasing function with radius. 

We assume a generic star formation history described by 
a gaussian function with a width given by the star formation 
timescale (tsp), which is left as a free parameter. The peak 
occurs at a time 2 x tsF and we assume an age of the Universe 
of 12 Gyr. Figure 3 shows the predicted relative abundances 
of this simple model as a function of tsF for both IMFs 
described above and for a set of ages: 0.5, 1, 5 and 11 Gyr. 
The number in parenthesis gives the slope of the gradient 
for both IMFs. It is beyond the scope of this paper to infer a 
proper functional dependence of the star formation timescale 
with radius. Hence, we will just assume a generic power-law 
dependence: 







(8) 



We find that the predicted abundance radial gradient for old 
stellar ages is A[Mg/Fe]/Ar ~ —0.15/3. In order to compare 
the model and the observations, we project the 3D radial co- 
ordinate (r) on to a 2D radius (R) using a generic luminosity 
profile: 

L(< r) oc r"^ ^ L(r) cx r^'^ (9) 
By projecting the volume distribution on to a surface, we 



can compare the observed slope a 
with the model prediction: 

a = 7- 2 -0.15/3. 



A [Mg/Fe] /A log i? 

(10) 

Using the deprojected function presented in Hernquist 
(1990), which gives a 2D de Vaucouleurs profile, we can 
choose 7 = 2. This value also matches a NFW profile 
(Navarro, Frenk & White 1997) for radii shorter than the 
core radius. In this case, the slope of the correlation between 
radius and star formation timescale can be written: 



/3(fixed 7) = -Q-^- 



(11) 



On the other hand, we can fix f3 much in the same way 
as most semi-analytic modellers do (e.g. Baugh et al. 
1998; Kauffmann, White & Guiderdoni 1993) so that the 
star formation efficiency (l/tsp) scales with the dynamical 
timescale ^rVGAf(< r). If we impose a mass-to-light ratio 
that does not change with radial distance — to be expected 
in galaxies formed in short duration bursts — we can use 
equation (^ to write M(< r) oc r~' , which implies: 

(12) 

and in turn implies: 
a + 2.225 

^ = ^:o7^- 

Our model assumes that the gradient of the abundance 
ratios is caused by a formation scenario with different burst 
lifetimes, which depend on the distance to the centre. One 
could argue that a large spread of efficiencies in different 
regions of elliptical galaxies would result in a rather large 
age spread. However, the age-sensitive Balmer absoption in- 
dex H/3 does not show any correlation with radial distance 
(Davies et al. 1993). This implies elliptical galaxies do not 
have a wide age spread in their stellar populations. Figure 4 
shows that even though we consider a large range of star for- 
mation timescales, our model is in agreement with the data. 
We plot in this figure the mass-weighted age of the stellar 
component as a function of galaxy age for three different 
star formation timescales. The inset gives the predicted ab- 
sorption in Hp as a function of stellar age for simple stellar 
populations with three different metallicities from the lat- 
est population synthesis models of Bruzual & Chariot (in 
preparation). One can see that after 3 — 4 Gyr, it is virtu- 
ally impossible to determine the stellar ages. This is the age 
shown as a dashed horizontal line in the figure, which implies 
that for galaxy ages of 8 Gyr, we can explore a range of star 
formation timescales up to tsp ~ 2 Gyr and still reproduce 
the observed lack of correlation between radius and Balmer 
index. 



5 DISCUSSION 

Our model links the abundance ratios to the star formation 
timescale, which cannot be directly observed. Instead, the 
data available give radial gradients. A simple power law, as 
described in equation (^), can be used in order to throw 
light on the possible correlation between the infall rate or 
the star formation efficiency and the dynamical properties 
of the galaxy. Figure 5 compares the observed radial gradi- 
ents with the predicted ones, as a function of [Mg/Fe]. The 
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6 8 10 -0.2 +0.0 +0,2 +0.1 +0.6 

Age/Gyr [Mg/Fe] 



Figure 4. Mass-weighted age of the stellar component as a func- 
tion of galaxy age for our generic star formation history. The inset 
shows the predicted equivalent width of Hf3 for several metallic- 
ities (from the top Z/Zq = {0.5,1, and 2}) using the latest pop- 
ulation synthesis models of Bruzual &: Chariot (in preparation). 
Mass-averaged ages older than (t)M ^ 3 Gyr (arrow) cannot be 
discriminated using Balmer indices. 



hollow squares represent the data from Trager et al. (2000a), 
measured at two projected radial positions: Re/8 and Re/2, 
which enables us to determine a gradient. A clear trend is 
seen towards positive gradients in more [Mg/Fe]-enhanced 
(i.e. more massive) galaxies. The other data points come 
from the compilation of Kobayashi & Arimoto (1999) from 
which we selected the work of Gonzalez (1993, GON); Davies 
et al. (1993, DSP), and CaroUo et al. (1993, CDB) estimated 
by the authors of the compilation to have the most reliable 
data. There is a very large scatter in this compilation, with 
a wide range of slopes, both negative and positive. In agree- 
ment with the data from Trager et al. (2000a), a correlation 
can be seen so that galaxies with positive slopes tend to have 
the largest super-solar abundances, which correspond to the 
most massive galaxies as seen in figure 2. The solid line gives 
a least-squares fit to the data from Trager et al. (2000a). 
The dashed line also shows a similar trend in the compila- 
tion of Kobayashi & Arimoto (1999). Estimating abundance 
ratios from line indices is still a very delicate issue espe- 
cially because of the difficulty in finding a suitable set of 
stars for calibration. Hence, a theoretical approach is called 
for, using the computations of Tripicco & Bell (1995), who 
recomputed all of the Lick/IDS spectral indices from a grid 
of theoretical spectra and atmospheres with varying abun- 
dance ratios. The response functions to nonsolar ratios is 
used to correct the standard population synthesis models of 
Worthey (1994), calibrated for solar abundance ratios. The 
uncertainties in such translation frmo abundances to line in- 



Figure 5. Observed abundance gradients estimated from the 
measurements of Trager et al. (2000a; TRG; hollow squares) at 
Re/S and Re/2. A typical error bar is shown for illustration. The 
solid line is a linear squares fit to this data. A rough estimate 
of [Mg/Fe] gradients from Kobayashi & Arimoto (1999) is also 
shown. The data points from this compilation correspond to the 
sample of Gonzalez (1993, GON); Davies et al. (1993, DSP); and 
CaroUo et al. (1993, CDB). The trend clearly seen in the data 
from Trager et al. (2000a) as well as hinted at by the void in 
the lower right corner of the data from Kobayashi & Arimoto 
(1999), and illustrated by the dashed line, can be interpreted as 
a hint towards a correlation between galaxy mass (which scales 
with [Mg/Fe]) and the lower limit of the gradient. 

dices is discussed in Trager et al. (2000a), and an error bar 
including these uncertainties is shown in figure 5. 

A negative gradient implies an inside-out formation pro- 
cess, which is suggestive of a monolithic collapse scenario, 
with gas falling towards the centres of dark matter halos, 
with a gas density profile that decreases outwards and a tem- 
perature that decreases inwards, thereby generating a pro- 
cess of star formation that starts at the centre and spreads 
outwards. On the other hand, recent TreeSPH numerical 
simulations of galaxy formation (Sommer-Larsen, Gotz & 
Portinari 2002) find early-type galaxies (both ellipticals and 
lenticulars) in which star formation proceeds in the oppo- 
site way, namely outside-in, mainly triggered by merging. 
Furthermore, recent observations of field spheroidals in the 
Hubble Deep Field (Menanteau et al. 2001) show blue cores, 
result which is currently interpreted as secondary bursts 
of star formation, but which could be an indication of an 
outside-in formation process, which implies a positive gra- 
dient. 

In the light of these data, we infer a varying correlation 
between the star formation timescale (tsp) and galaxy ra- 
dius. A simple approach, fixing 7 as in equation ( |ll| ) gives 
a slope for this correlation in the range —1.3 < fi < 2. A 
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better approximation, relating /3 and 7 through the dynami- 
cal timescale, assuming a fixed mass-to-light ratio (113) gives 
0.4 < P < 0.6 with a mass or luminosity profile slope of 
1.8 < 7 < 2.2. The higher values of 7 correspond to positive 
gradients in the abundance ratio. 

The scatter is rather large, and it can only imply a 
very weak correlation between /3 or 7 and some dynami- 
cal parameter such as central velocity dispersion or mass. 
This scatter gives another hint of merging as the major 
mechanism in the assembly process of early-type systems. 
In a purely monolithic collapse scenario, one would expect 
a well-defined correlation between the dynamical timescale 
and the star formation timescale, in such a way that most of 
the early star formation would occur at the centre. The void 
in the lower right corner of figure 5 (illustrated by the thick 
dashed line) shows that we can infer a correlation between 
[Mg/Fe] (or galaxy mass as shown in Figure 2) and a lower 
limit to the gradient. This can be interpreted as monolithic 
collapse being possible only in low mass systems. Higher 
mass galaxies must be assembled through merging. We em- 
phasize here that the inside-out versus outside-in classifica- 
tion may be an oversimplification. Depending on the simple 
model we choose — i.e. fixed luminosity profile, as in (|ll|), or 
correlating /3 with the dynamical timescale, as in (|l2|), mas- 
sive galaxies give respectively real outside-in star formation 
(/9 < 0) or a shallower slope for the correlation between star 
formation timescale and radius. In either case, the correla- 
tion suggested by the dashed line in figure 5 hints at more 
star formation in the outskirts of massive galaxies compared 
to low mass systems during the major stage of star forma- 
tion. 

It is worth mentioning that the slope of the abundance 
ratios is a more robust estimator than the absolute abun- 
dance. The latter is strongly dependent on the amount of 
gas ejected in outfiows, which is an important mechanism in 
early-type galaxies as hinted at by the mass-metallicity rela- 
tion (Arimoto & Yoshn 1987; Ferreras & Silk 2001), whereas 
the slope would depend on the scaling of the outflows with 
radial distance, which should be similar to the scaling of 
the star formation efficiency and infall timescale, thereby 
reducing the effect of outflows on the slope. 

The radial gradients in abundance ratios represent an 
alternative way of studying the universality of the IMF, com- 
pared to analyses of the metallicities of local stars or the in- 
tracluster medium ( Wyse 1997) . Figure 3 shows that a non- 
universal IMF would translate into a slope change of the 
abundance ratio [Mg/Fe]. For instance, if we expect a top- 
heavy IMF in environments with a high star formation rate 
(i.e. with a short tsp) then the correlation between [Mg/Fe] 
and star formation timescale should be steeper. However, 
this calls for a more detailed model which is beyond the 
scope of this paper. We have also neglected the effect of 
early galactic winds which could eject a-enhanced material 
out of the galaxy. A correlation of these winds with the local 
escape velocity could explain the observed radial gradient of 
Mg abundance in ellipticals (Martinelli, Matteucci & Co- 
lafrancesco 1998). 

The model presented here relies on the fact that we un- 
derstand the mechanisms that trigger SNIa and we are able 
to predict their rates or at least their ratio with respect to 
SNII. Short bursts of star formation imply that most of the 
contamination from SNIa will go to the ISM and not so much 



to the stellar component. Hence, we expect gas in ellipticals 
to have solar abundance ratios. There has been a long con- 
troversy over this point. The study of ASCA observations 
of a few clusters hinted at a signiflcant SNIa contamination 
of the gas in rich clusters (Ishimaru & Arimoto 1997). How- 
ever, Gibson et al. (1997) showed that a different calibra- 
tion could invalidate this hypothesis. Arimoto et al. (1997) 
even raised the doubt of whether abundance estimates using 
the standard iron L-line complex in X-rays is giving wrong 
metallicities. The latest analysis of six early-type galaxies in 
Virgo using ROSAT and ASCA observations (Finoguenov & 
Jones 2000) seem to agree with a high iron content, i.e. an 
important contribution from SNIa to the enrichment of the 
ISM in elliptical galaxies. 

In this paper, we conclude that a simple treatment of 
the abundance ratios allows us to infer the star formation 
history and its connection to the dynamical formation his- 
tory. The data seems to rule out a monolithic or secular 
formation scenario for the stellar component of massive el- 
lipticals. On the other hand, the large scatter in the slope of 
the radial dependence of [Mg/Fe] in low-mass galaxies im- 
plies both a merging and a monolithic "mechanism" should 
be invoked for these systems. A combined analysis of abun- 
dance ratios both in the stellar populations and in the in- 
terstellar medium will enable us to explore the connection 
between these two histories, and to quantify the importance 
of merging events both in the star formation and dynamical 
histories of early-type galaxies. 
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